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Abstract 



We investigate statistical properties of the eigenfunctions of the Schrodinger op- 
erator on famihes of star graphs with incommensurate bond lengths. We show that 
^ ' these eigenfunctions are not quantum ergodic in the limit as the number of bonds 

tends to infinity by finding an observable for which the quantum matrix elements do 
not converge to the classical average. We further show that for a given fixed graph 
QQ . there are subsequences of eigenfunctions which localise on pairs of bonds. We describe 

' how to construct such subsequences explicitly. These constructions are analogous to 

scars on short unstable periodic orbits. 

o 

^! 1 Introduction 

J2 ' Let ^l^n denote the wave-function corresponding to the n^^ energy level of a quantum 

system that has a Hamiltonian dynamical system as its classical limit. We are interested 
^ ! in these wave-functions in the n ^ oo limit, which corresponds to the semi-classical regime. 

I Numerical and some theoretical evidence supports the hypothesis that their behaviour in 

^ ■ this limit is determined by general properties of the underlying Hamiltonian such as, for 



example, time-reversibility, integrability and statistical properties of the flow (ergodicity, 
mixing, etc.). A deeper understanding of this is one of the goals of current research in 
quantum chaology. 

When the classical Hamiltonian generates chaotic motion, the semi-classical eigenfunc- 
tion hypothesis asserts that the wave-functions should equidistribute over the appropriate 
energy shell |BeH fV] . A physical explanation for this is that in the semi-classical limit the 
quantum system should mimic the behaviour of the classical system; if the classical motion 
is chaotic, then a typical trajectory ergodically explores the surface of constant energy in 
phase space. Another interpretation is that eigenstates are invariant under time evolution, 
so it is natural to associate them in the semi-classical limit with classical invariant sets. 
One such invariant set is the energy shell itself. 

The Schnirelman theorem ICdV| [Zj KirLj states that for systems in which the Hamil- 
tonian flow is ergodic the sequence of measures induced by ipn converges to Liouville 
measure in the limit as n — > oo along a subsequence of density one. This behaviour has 
been termed "quantum ergodicity". Quantum ergodicity implies a weak version of the 
semi-classical eigenfunction hypothesis ^ 



It is possible that quantum ergodic systems have subsequences of states for which the 
corresponding measures do not converge to Liouville measure (of course such subsequences 
have density zero). These subsequences, if they exist, are expected to be associated with 
other classical invariant sets, such as periodic orbits. The case where the limit of an 
exceptional subseqence is a singular measure supported on one-or-more isolated, unstable 
periodic orbits of the classical system is called "scarring" . 

Scarred eigenstates were observed numerically by Heller [H|, who proposed the first 
theoretical explanation for their existence, based on the semi-classical evolution of a wave- 
packet centred on a periodic orbit under linearised dynamics. Another important develop- 
ment was an understanding of the contribution to wave-functions from all periodic orbits 
|Bol IBe2| resulting in formute related to the semiclassical trace formula for the density 
of states. Later, the theory was extended to include non-linear effects |KHj and, more 
recently, situations where the orbit in question undergoes a bifurcation |KPj . A review of 
related works was given in |Klj . All of the above mentioned theories relate to scar effects 
in averages over a semiclassically increasing number of states. This may be thought of as 
a weakened form of scarring, because it is not clear that any one state in the averaging 
range causes the scar; the scars may be a collective effect. It is a much harder problem to 
show that a particular sequence of individual states is scarred. Currently, the only systems 
known rigorously to support scarring in this strong form are the cat maps |FNdB] which 
have non-generic spectral statistics caused by number-theoretical symmetries |Kej . 

For quantum graphs |KSlj the wave- functions are the eigenfunctions of the (continuous) 
Laplace operator on the bonds with matching conditions at the vertices chosen to make 
the problem self-adjoint. There is evidence to suggest that the spectral statistics of large 
quantum graphs coincide with those of generic quantised, classically chaotic, systems |KS2l 
IBSWHfBSW2| lB] subject to mild conditions on the connectivity T . Although the classical 
dynamics on a graph are not Hamiltonian, they are ergodic and it might be expected that 
an analogue of the Schnirelman theorem should hold. 

Graphs are a rich source of problems in quantum chaology and related fields. Re- 
cent works have considered: scattering problems .KS3> .TM, KS4 , the spacing distribu- 
tion of eigenvalues |B(Tj . nodal domain statistics [(jSWj . the Dirac operator on graphs 
|BH11 IBH2j , Brownian motion on graphs |CDM1 ID] and the important question of how to 
construct families of graphs with increasing numbers of bonds jPTZj. 

Recently, authors have begun to investigate the wave-functions of quantum graphs. 
Kaplan |K2j studied eigenfunction statistics for ring-graphs using a combination of nu- 
merical techniques and analytical calculations of the short-time semiclassical behaviour of 
a wave-packet close to a 1-bond periodic orbit. The inverse participation ratio (a mea- 
sure of localisation in a given state) was found to be well-described by this contribution, 
and shows deviation from the ergodically expected behaviour. Similar deviations were 
noticed for lattice-graphs. Remarkably, Schanz and Kottos |SKj observed that it would 
be impossible for the shortest orbits that are responsible for this enhanced localisation to 
support strong scarring. They wrote down an explicit criterion which must be satisfied 
by the energy of any strongly scarred state, and deduced asymptotics of the probability 
distribution of scarring strengths. In |KMWj a study was made of the eigenfunctions of a 
family of graphs known as star graphs (the name being derived from the connectivity of 
graphs in the family). The value distribution for the amplitude of eigenfunctions on a sin- 
gle bond of the graph, subject to an appropriate normalisation, was rigorously calculated 
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in the limit as the number of bonds tends to infinity. In fact the normaUsation impUes that 
star graphs with a fixed, finite number of bonds are not quantum ergodic. However, this 
result leaves open the question of whether star graphs are quantum ergodic in the limit as 
the number of bonds tends to infinity. This is because one bond represents a vanishingly 
small fraction of a graph when the number of bonds becomes infinite, whereas quantum 
ergodicity is concerned with structures on macroscopic (classical) scales. 

The results we present here extend the work in |KMWj on star graphs. We review the 
definition of a quantum star graph in section |2l below. We show that (see the following 
subsection for precise statements) quantum star graphs are not quantum ergodic in the 
limit as the number of bonds tends to infinity. We also show that for any given star graph 
there exist exceptional subsequences of eigenfunctions that become localised on pairs of 
bonds as n ^ cxD. Orbits on a graph are simply itineraries of bonds, so this localisation 
is analogous to strong scarring on short period-2 orbits. Such orbits are unstable in the 
sense that there is an exponentially small probability of remaining on a given orbit. Our 
explicit construction supports the observation of Schanz and Kottos [ SKI that star graphs 
support a large number of states scarred in such a way. 

The spectral statistics of star graphs are different to those associated with the more 
general graphs described above |BKj . The fact that quantum star graphs are not quantum 
ergodic does not contradict the possibility of a quantum ergodicity theorem for graphs 
with general connectivity. It is known that the spectral statistics of quantum star graphs 
are the same as those associated with the family of Seba billiards |Se[ IBBKj , so-called 
"intermediate statistics". There is evidence to suggest that the results we present on 
scarring can also be extended to Seba billiards |BKWj . 



1.1 Main results 

To investigate quantum ergodicity for large star graphs, we consider an observable that 
picks out a positive proportion of the graph. We consider a graph with av bonds, where 
a,u G N, and the observable B = {Bi{x))f^-^ defined by 

Br.= \l = (1.1) 

[ (J tor I = V + 1, . . . ,av. 

B may be thought of as the indicator function of the first v bonds. The classical average 
of B is approximately 1/a. We shall consider the limit v oo. 

Wave-functions on graphs have a component on each bond, so we shall use the notation 

lA^'^) := (#^)r=^ 

for the n}^ eigenstate. The inner product (•!•) is defined in (|2.8|) below. 

Each bond of the graph has a length, and the vector of bond lengths will be denoted 

Theorem 1.1. For each v let the components ofL he linearly independent over Q. Then 
there exists a probability density Pviv) such that for any continuous function h, 

hm J^/i((V,(")|B|tA(")))= / /i(??)p.(??)d??. (1.2) 

The density py{ri) is supported on the interval [0,1]. 
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Theorem 1.2. For each v let the bond lengths Lj, j = 1, . . . ,av lie in the range [L, L+AL] 
and be linearly independent over Q. If vAL as v ^ oo then there exists a probability 
distribution function F{R) such that for any R £ (0, 1), 



and 



lim 

I)— >oo 



R 



(1.3) 



where 



/ TTu; J — oo 



ri=l/R-l 



(1.4) 



TTT] 



exp 



-ivr (a-1) /ivr i^^ 



+ ^^^^exp(-^-^) +e(a-l)erf 



The function F{R) is plotted in figure E 

Remark 1.3. // star graphs satisfied quantum ergodicity, then F{R) would be the step- 
function 

^(^) = l n {"^^'^w" (1-5) 
^ ^ { 0, for R < 1/a ^ ' 

for this observable. 

In figure n we compare the numerical data for the value distribution of (^''"'•*|B|^^"^) 
for a star graph with 90 bonds with the v ^ oo analytical prediction F{R). The difference 
between the actual distribution F{R) and that which would be expected if the graph were 
quantum ergodic (remark \l.'6\i is clear. Figure [21 shows the difference between numerical 
data and F(R) for increasing values of v. 

We we also show that for graphs with fixed number of bonds, there are subsequences 
of eigenfunctions that localise on two bonds. 

Theorem 1.4. Let the elements of Jj be linearly independent over Q. Given any distinct 
two bonds, indexed by ii and i2, of a v-bond star graph, there exists a subsequence (kn^) C 
(kn) such that for any f = {fi)i^i smooth in each component, 

lim (^("^^IflV^"'-)) = , ( ^'^ /n(x)dx + f^,{x)dx) . (1.6) 
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Figure 1: Comparing numerical data with the analytical prediction, F{R). For this plot 
a = 3 and in the numerical study, v = 30. 
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Figure 2: Convergence to for v = 5(+), 10(x), 15(*), 20(Q), 25(H), 30(0). 

2 Quantum Star Graphs 

A star graph^ is a metric graph with b vertices all connected only to one central vertex. 
Thus there are b + 1 vertices and b bonds (figure OJ- We shall denote by L G M'' the vector 
of bond lengths. 

We define the quantum star graph in the following way. Let TC denote the real Hilbert 
space 

n:=L\[0,Li]) X... xL\[0,Lh]) (2.7) 

with inner product 

(f|g):= fj{x)g,{x)dx. (2.8) 

^sometimes referred to as a Hydra graph 
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Figure 3: A star graph with 5 bonds 



Elements of 7i are denoted f = (/i, . . . , /{,). Let J- ^Tihe the subset of functions f which 
are twice-differentiable in each component and satisfy the conditions 

/i(0)=/i(0) =: /o, j,i = l,...,6 (2.9) 

^ 1 

f,{Lj) = 0, j = l,...,6. (2.11) 

The parameter A may be varied to give different boundary conditions at the central vertex 
of the graph. Henceforth we shall concentrate on the case 1/A = 0, the so-called Neumann 
condition. The Laplace operator A on .7^ is defined by 

A defined on T is self-adjoint. Since the space on which the functions in are defined 
is compact, the operator A has a discrete spectrum of eigenvalues f |DSj . Section XIII. 4). 
i.e. the equation 

-AV' = /c^i/; (2.13) 
has non-trivial solutions for k = ki,k2,--- Such ^ are the wave- functions |KSll IK^ . 
We shall use the notation that := {ipl"'\x))^^i is the wave-function corresponding to 
k = kji- 

Solving 1)2. 13(1 with boundary conditions 1)2. 9(1 - 1)2. 11|) . we find that the component of 
the n*^ normalised eigenfunction of the Laplace operator on the z*^ bond of a star graph 
is 

V'f ^(x) = cos kn{x - Li) (2.14) 
where the amplitude is given by 

and kn is the n^^ positive solution to 

b 

Z{k,L) ■.= ^tankLj = 0. (2.16) 
i=i 

In sections OHSl it will be convenient to take b = av, where a € N is fixed. This is so 
that we can easily describe a fraction of the total number of bonds as the number of bonds 
becomes large {v — > oo). In section |H1 we shall take a = 1 for notational convenience, since 
there we will only be concerned with fixed graphs. 
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3 Distribution of the observable B 



In this section we prove the existence of a hmit distribution for the diagonal matrix 
elements of B on star graphs with a fixed number, av of bonds. 

Lemma 3.1. Consider a star graph with av bonds with fixed lengths given by the vector 
L. Then for B defined by (jl.lj) . 

iWlRlW.Wv - Ei=i sec2 /c„Li , ^ 1 



Thus 



where the error estimate is uniform in v and Li ^ Lmin > for each i. 
Proof. We recall that 

av 

i=i 

Integrating ()2.14l) gives, for i = 1, . . . , f , 

/ \i{jI"\x)\'^ Bi{x)dx = — 21, T {LiSec^KLi + ^tanknLi], (3.19) 
and for i ^ V + 1, 

Li 

|^f^(x)|2Bi(x)da; = 0. (3.20) 

= + ^ (3.21) 

2l^j=X -t^j sec Kn^j n^ra 

where 

X]i=i tan fcnLi 
~ V"" L sec2fc L ' ^ ^ 

Let Ljnin miuj {Lj}. Then 

l^gij ^ Ya=i I tanfcnLtl (3 23) 

(a - l)Linin^^ + -i^min I]Ll SeC^ 

< (l...,„+ '"-"fr; )" (3.24) 

using the fact that | tan0| ^ sec^ 6 for any ^ G M. Hence E = 0(1) as n — > 00 uniformly 
in V, Lmin > 0. □ 

Proof of theorem By lemma 13.11 

M(V(")|B|i/.("))) = h ( ^b^/''"y ] + En (3.25) 

where En = o(l) asn— >oo since h is uniformly continuous on [0, 1]. Hence 

1 ^ 

— ^ £;„ ^ as ^ 00. (3.26) 

n=l 
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Therefore 



lim lf;M(V(-)|B|V'('^))) = lim if;/. ■ (3.27) 

n=l n=l — ■> ■> / 



1 ^ r 
J™ l^E/(^"L)= / /(Odz^(0 

n=l ^ 



According to Barra and Gaspard, there is an absolutely continuous measure vi^ such 
that for piecewise continuous functions, / : S ^ M, 

i ^ 
m — 

N- 

where S is the surface embedded in the av dimensional torus with side vr, defined by 

tan xi + • • • + tan x^v = 0. 

^ is a set of af — 1 coordinates which parameterise S. To avoid repetition, we refer the 
reader to |BGj . |KMWj for more detail about this result and its application to similar 
problems. 

Let 

/(x) = h f Siilf^^ (3.28) 

we can define p„ {r\) by 

/oo 
/i(r?)p,(r/)dr?. (3.29) 
-oo 

Since ^ (■0(")|B|i/)(")) ^ 1, it follows that p„(7?) is supported on [0, 1]. □ 

4 The large graph Umit 

Let ry G M and define 

^ CIV V 

Xrj{n) := — ^ Lj sec^ A;„Lj - -^'^Li sec^ A;„Lj (4.30) 
for n = 1, 2, . . . The key result of this section is the following. 

Proposition 4.1. For each v, there exists a probability density function fx,^,v such that 
1 

hm # {n G {1, . . . , AT} : X,(n) < 5} = / /x„.(tT)da. (4.31) 
Furthermore, for each S G M, 

/ /x^,^(o-)dCT ^ / /x^(f7)d<T 
as V ^ oo, provided that vAL in this limit, where 

The functions Pr^ andrrj are defined by (|4.43|) and (|4.44j) below, andw{z) := e~^^ erfc(— iz). 
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Proof. The existence of the N ^ oo limiting density, fxn,v is a consequence of the result 
of Barra and Gaspard |BGj . The proof is entirely analogous to the proof of theorem 11.11 

We turn our attention to the v ^ oo limit of this density. With n a random variable 
uniformly distributed on the set {1, . . . , A^} for some N £ N, the characteristic function 
for this random variable Xn{n) is 

1 ^ 



n=l 



where / : [0, vr]^ ^ C is defined to be 



— sec^ Xj — rj sec^ Xi 

Following the argument of |KMWj we can write 
1 ^ 



n=l 



1 1 />00 PTT PIT I \ I -A \ 

j J - J I ^^^^ 1 •^(^) I 7 1 d^^'xdC. 



Denoting this limit for each fixed by e^(/3), we can write 

e.(/?) = ^ y /i/r'^s"-" + (« - l)/4/2"/3"'^-"-MC (4.32) 

where the integrals /i, . . . , I4 are:- 

/,:=irsec=.exp(l«.a„.-l^.ec=.)d., (433) 

I2 '■= — [ exp ( — tanx — sec^ X ) dx, (4.34) 

An integral similar to (|4.32|) was tackled in |KMWj . We quote here the relevant results, 
mutatis mutandis. Asymptotic analysis of the integrals in (|4.33M1^!^ gives 

1' / ivr iC"^ \ 

exp -— + -|- +0(1) asi;^oo (4.37) 



and 



V /ivr iC^\ ^ / s / ^ 

^exp -t-t;? +0(1) asi;^oo. (4.38) 



vr/? V 4 4/3 
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For I2 and /3 we consider separately the cases —y/v < C, < y/v and \C,\ > \/v. For 
and 

4"-"— (-^^-(-T-S) ©)) (-) 



both error estimates are as w — 00. For |C| > aA', 



as ^ — ^ 00, and 

|/3K^f4 + ?)'' + 0(rt (4-42) 

Using these estimates, we can find an expression for the hmit of e„(/3) as « ^ 00, 

40) := mn e„(« = ± £ -Lp, (£) exp (-^^r, (^)) dC 

Pv(()<'^p(-/PT,(())d(. 

For ease of notation, we have introduced 

and 



r.(0 :~V^exp(^ + |-)+eerf 



ivr , i^2^^ , . cf 



+ exp (- J - t J + - 1) erf j . (4.44) 

In the above, wherever occurs for /3 < 0, this should be understood to mean ±i-\/— 
where the sign is taken in such a way that 



e{-(3) = e{f3), 

the usual condition for the characteristic ftmction of a probability density. This can always 
be done. We also note that e(0) = 1 which is consistent with e{0) being the characteristic 
function of a probability distribution. e(/3) is continuous at /3 = since the defining 
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integral is uniformly convergent in (3 (see lemma lA . 31 below) . Thus the limiting density, 
exists and is given by 

/^.(^) = ^/ ^ P,(Oexp(-v^r,(e)-ia/3)ded/3, (4.45) 

where we have made the substitution ^ = C/\^- We here switch the order of integration. 
This is a non-trivial operation since both integrals are improper. However in this case we 
can rigorously justify the manoeuvre. Justification is provided in appendix El proposition 

roi 

I f'OO f'OO 

fx,{<y) = / PviO exp(-v^r,(0-icT/?)d/3de 

- 1 r p,,o ( > _ 4f^. { s!^) (4.46) 



2-Ka ^ 2 (-a)3/2 



To evaluate the final integration we have used the following result (a variant of formula 
3.462.5 in [HEl), 

/ exp(-ax-6V^)dx = i-^-|^^(^). (4.47) 
To conclude we observe that since 

fOO 



/oo 
-oo 



the first term in 1)4. 46() vanishes as it has zero real part. □ 

Some properties of w{z) are discussed in |ASj chapter 7. We shall use the following 
Lemma 4.2. The function w{z) has the asymptotic expansion 

, N i (2m,)! /, , X 

~ 7; E -^^^ (4-48) 

^ m=0 

as z ^ oo, valid for —— < arg z < —-. 
! ■'4 4 

Proof. This follows from the asymptotic expansion of erfc, 

~" (2m)! 



^fnz^^ erfc(2;) ~ 



{4z^)"'m^-' 



m=0 

as z ^ oo, I argz| < 37r/4, taken from AS^ (formula 7.1.23; see also iBlHal exercise 3.11). 
The series for w comes from making the substitution z —iz. □ 

Since erfc and w are analytic functions, they are bounded in the domain of validity of 
their asymptotic expansions quoted in lemma ^21 

Proposition 4.3. Let z G C and — 7r/4 < argz < 57r/4. Then 

l" M^PHP = f - =|> + ^ log N + iC« + O (^) (4.49, 

where Cr G M is independent of z, hut may depend on R. 
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z\R 



Figure 4: The contours 71, 72 and 'Jr. 



Proof. Write z = |z|e"'' where (p = avgz. Then 



R i-\z\R 

zw{zp)dp= / e''^w{e''^p)dp (4.50) 







via p p/\z\. Using Cauchy's theorem, 

/' ' e''^w{e''^p)dp = I w{t)dt = j w{t)dt + / w{t)dt. (4.51) 

The contours 71 , 72 and in the complex t-plane are illustrated in figure 0] 
On 72, 

w{t)dt = / w{x)dx 



72 Jo 



1 MR / ■ X i-\z\R J 

w{x)dx + / I w{x) ^=r- ] dx + dx 

Jl \ VTTXy J I ^TXX 

1 ,00 1 ■ \ 

'w{x^dx + / w{x) -j=^ dx 







00 , J 



^(a;) ) dx + / dx 

°° _ 2 /"°° / i \ /■''''^ i 

e ^ dx + i X const — / w{x) dx + / dx, 

J\z\R V V^2;y 7i ^-Kx 

using tt;(2;) = e~^ (1 + ierfiz) to separate the real and imaginary contributions. The 
imaginary part goes into the constant C/j, the value of which is not important since we 
shall always be considering only the real part of resulting expressions. By the use of lemma 

y ti;(t)dt = ^ + ix const + -J=(logi2 + log|z|) + o(^^^^ (4.52) 
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as R ^ oo uniformly for \z\ ^ c for some c. 
On 7ij, 

f(t)dt = r i\z\Re'''w{\z\Re''')d9 

IR Jo 



"^ilzlfe'^ ( ]^ n + o( . .1^.. ] ] bv lemma Ol 



y^/^i\z\Re'^ \\z[^R^ 



Combining and gives □ 

The following lemma from probability theory will also be useful. 

Lemma 4.4. Let U, V be random variables, then 

where fx^ is the probability density function of the random variable := U — r]V. 
Proof. This follows immediately from the fact that 

^ <r]] = F{U-r]V < 0). 



V 



5 Proof of theorem 11.21 

We first observe that 



□ 



Yli=i Lj sec^ knLj _ [ , Y^'j=v+i 



Lj sec k^Lj \ 



so we can concentrate on finding the probability distribution, F[rf), of 



(5.54) 



Xyl = l -^i S6C^ k^Li 

as f ^ cxD. This will be then related to the distribution in which we are interested by a 
simple transformation (see equation 1)5. 57() below). In light of lemma 14.41 F{r]) is given by 

fO 



J — oo 



with fx^ found in proposition 14.11 It is more instructive to take the range of integration 
from —R^ to 0, with a view to taking the limit i? — > cxd later, 



/x,(cT)da=^iHe r r ^.(0^^^^ f ^^^^^V^d^- (^-^^^ 
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The ,^-integral is uniformly convergent by lemma IB.ll so we have legitimately switched 
the order of integration. We can then write 



\/vra J-ooJ-oo 4(-ct)3/2 \^ j 

The second term vanishes as i? ^ oo, as shown in DroDosition lB.2l Making the substitution 
= 2p, leads to 

F{n) = / / P,(e)e3-/V,(0ti^(e=^"'/V,(e)p)dpde 

We can apply proposition 14.31 with z = e^"^/^r^(^). We can integrate ^ out of the error 
term provided by this proposition since 



d^ < OO. 



We get, finally, 



= i + ±9^e (arg(T,(e)) - ilog |T^(e)|) d^. (5.56) 

This is related to F{r]) by 

F(77) = 1 - F f 



F{r]) = 1 - F ( i - 1 ) . (5.57) 



□ 



6 Scarred states on finite star graphs 

We recall that theorem 11.41 is concerned with the quantity (t/)^"^|f |^^"^) which can be 
written as 

^ rL. 



Now using (I^TH) and the identity cos^ 6^ = 1/2 + (1/2) cos 26 yields 

(^W|f|^W) = ^ -1^ i^J^ fi{x)dx + cos2A:„(x - L,)/,(x)dxJ . (6.58) 



Since we are interested in a subsequence — > oo, we may hope that the second integrals 
do not survive. Thus our prime concern are the prefactors 



ELi Lj sec2 knLj ' 
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and to prove the theorem we need to find kn such that the prefactors are small for all i 
other than (without loss of generality) 1 and 2. 

To do this we study the poles of the function Z(A;,L). Let 

vr /I 



Then | tan pLi\ oo and sec pLi — > oo as p — > pn^i- 

Since the function Z{k, L) is an increasing function of k, between any two poles there 
is a zero. We will use this important feature to "trap" zeros k^^ of Z{k, L) between pairs 
of nearby poles pn^i and Pm,2, also requiring that all other poles are far away (see Fig. [nj. 
The implications will be that as r — > oo, 



and 



sec^ kn^Li ^ sec^ kn^Lj with i = 1,2 and j > 2, 



sec^ kn^Li ~ sec^ kn^L2, 



ensuring that (|1.6() holds. 

To make the above arguments rigorous we need the following propositions. 

Proposition 6.1. Let the elements of Jj be linearly independent over Q. Let 1 < v* < v 

for V ^ 3. Given e > there exist infinitely many n G N such that 

• for each i = 2, . . . ,v* there exists m G N satisfying 

\pm,i -Pn,i\ ^ e/2 (6.59) 

• for all i = V* + 1, . . . ,v and for a// m E N 

vr 

\Pm,i - Pn,l\ ^ ^/^- {6.60) 

Proof. The idea behind the proof is that for linearly independent elements of L the poles 
p.^i for different i behave like independent random variables, therefore every permitted 
pole configuration happens infinitely often. To substantiate this claim we express the 
nearest-pole distances as the states of an ergodic dynamical system. 

For n G N and i = 2, . . . ,v* , let 5n,i denote the distance between and the closest 
pole of tan kLf, 

^n,i . — Pn,l Pm,i) 

where m is such that 

\Pm,i -Pn,l\ = vahl{\pm,i -Pn,l|}- 
m 

Since the poles of tan/cLj are vr/Lj-periodic, we have 

TT TT TT 

^n,i + Trr = Pn,i-Po,i + -^ mod— (6.61) 

TT TT TT 

H n mod — . (6.62) 



2Li Li Li 
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r 



Figure 5: Poles p^^j and nodes on the real line. Different symbols correspond to 
different values of i (the circle corresponds to Pn,i)- Filled symbols correspond to the 
poles, empty symbols to the nodes. In this example u = 5 and = 2. 



Let im,i denote the vn!'^ zero of tankLi. We note that (|6.60j) is implied by the condition 
that 

for some m G N. For i = v* + 1, . . . ,v, define rjn^i to be the distance between pn,i and the 
closest zero of tankLi. Similarly to ()6.62|) 

vrvr/llXvr vr 



From (|(1H2|) and 6n,i and r]n,i satisfy the recurrence 

6n+i,i = 6n,i + mod f- i = 2,...,v* 
r]n+i,i = Vn,i + fi mod f- i = v* + l,...,v 



(6.64) 



Since the bond lengths are not rationally related, the dynamical system ()6.64() is equivalent 
to an irrational translation on a torus. In this case, Weyl's equidistribution result |Wj 
applies, and any subset of the torus with positive Lebesgue measure is visited infinitely 
many times. The volume of the area in 5 — r/ space defined by 

-e/2 <5n,uVn,i<e/2 (6.65) 

is non-zero and so there are infinitely many n for which ()6.65() . and therefore 1)6. 5916. 6U() . 
are satisfied. □ 

The interpretation of proposition 16.11 is that we can find situations on the real line 
where v* poles of the functions tankLi are bunched together and the remaining v — v* 
poles are not close to these bunched poles (see figure 1^1). 

Proposition 6.2. Under the conditions of vrovosition f^TTl there is a subsequence (kn^) C 
(kn) for which 

sec^ kn^Li — > oo for i = 1, . . . ,v* 
sec^ kn^Li — > 1 for i = v* + 1, . . . ,v 

as r ^ oo. 
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Proof. Let (e.^) be a sequence satisfying e.^ ^ as r ^ cxd. We choose fe„^ as follows. 

Applying proposition IH. II with e = Er yields a set of v* poles of Z{k, L) inside a region 
with width e^- Since there is a zero of Z(A:,L) between any two poles of Z{k,'L), we can 
find i;* — 1 zeros in this region. Set A;„^ to be one of these zeros. 

From proposition we have 

\knr ~ Pm,i\ ^ Er for alH = 1, . . . , f * and some ra = m(r, i) 
\knr — ^m,i\ ^ Er for alH = i;* + 1 , . . . , u and some m = m{r, i). 

Since sec^ LjPm.i = OO; sec"^ Lii^^i = 1 and sec9 is a periodic function, the statement of 
the proposition follows trivially. □ 

Corollary 6.3. If v* = 2 in vrovosition \6.'A then we additionally have 

sec^fcn^Li 

^™ 2 U T = 

r^oo sec^ Kn^L2 

Proof. We recall that since kn^ is an eigenvalue, Z{kn^,'L) = 0, and hence 

tankn^Li = —ianknrL2 — tanA;„^L3 — • • • — tankurLv- (6.66) 



On the other hand, by proposition 16.21 tanA;„,,Lj remains bounded for i > 2 and tends to 
infinity for i = 1,2. Dividing (|6.66j) through by tan kn^L2 we obtain 

tanA;„,,Li 
lim ; — — = —1. 

r-+oo tan Kn,.^2 

Further observations that sin^ kn^Li — ^ 1 for i = 1, 2 and sec^ 6 = tan^ 0/sin^ conclude 
the proof. □ 

Lemma 6.4. Let / : [0, L] — > R be continuously differentiahle. Then 

f-L 



k—foo _ 

Proof. Integration by parts yields 



lim / cos(A;x)/(x)d3; = 0. 
Jo 



/ cos(A;x)/(x)d2; = — ( sin(kL)f(L) — [ sm(kx)f'(x)dx 
Jo k \ Jq 

and the statement follows immediately from the boundedness of / and its derivative. □ 

Proof of theorem \1.4i Without loss of generality, we can assume that ii = 1 and 12 = 2. 
We take the subsequence whose existence is guaranteed by proposition 16.21 with v* = 2. 
By corollary 16. 3| 



lim A 



M^-J^r. sec2V^\ _ f 2{L, + ii i = 1,2 




r^oo sec^ kn,.Li j [ otherwise. 



We use lemma ins to get rid of the second integrals in 1)6. 58(1 and conclude 
lim (V'("')|f|r/'("'-)) = ( h{x)dx + f2ix)dx) . 



17 



□ 



Acknowledgements: GB is grateful to the University of Bristol for hospitality during 
visits while part of this research was carried out. BW wishes to thank the University of 
Strathclyde for hospitality. 

BW has been financially supported by an EPSRC studentship (Award Number 0080052X). 

We gratefully acknowledge the support of the European Commission under the Re- 
search Training Network (Mathematical Aspects of Quantum Chaos) HPRN-CT-2000- 
00103 of the IHP Programme. 



A Appendix: The order of integration in (14.451) 

In this appendix we deal with some technical issues regarding the exchange of order of 
integration in (|4.45|1 . 

We first consider some asymptotics of Trj{S,). 
Lemma A.l. For ^ G M, 

r,(6=«ICI + 0,(r') 
as 1^1 — > oo, where the error estimate depends on r/. 

Proof. We first note that is an even function, so we may assume ^ > 0, and the result 
for < will follow by symmetry. We can write t^(^) as 



T,(£) = Vvt ( +{a- rm 



where 



We expand the error function asymptotically, 



erf = 1 — eric 



as ^ — > oo. Substituting HA.68|) into ()A.67|) gives 

t(e) =C + 0(r'), ase^oo, (A.69) 

and the lemma follows. □ 

dr 

Lemma A. 2. For ^ > 0, 9^e ^ and for all ^ G M, there exists r* > such that 
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Proof. By differentiation, 



di 



erf 



2 ; 



i7r/4 



vr 



5/2 



e""^ dr. 



We see that 



Thus, 



9le 



dt 2 /-^/^ 



cos r 



vr 



dr > 0. 



dT„ 



+ (a - l)9^e ^ 0. 



Hence it follows that fKe r^(0 ^ «Ke r^(0) = \/2(^ + a - 1)/V^ =: r* 
Lemma A. 3. The integral 



□ 



P,(0exp(-V/3r^(0)de 
is uniformly convergent for (3 G [0, /3o] /or a// /?o > 0. 

Proof. By making the substitution = we can consider the uniform convergence of 



Pr^i^/v) exp(-v^r^(Vi^)) 



di^ 



Let 



/(i/,/3) := J^exp(-v^9ler^(^)), 
0(1., /?) := :^^exp(-v^iamr,(^)). 



By lemma IX!T1 Jm t^(\/^') = 0(z. ^) as — > oo (We drop the r/-dependence since we are 
concerned here only with fixed rj). So 



exp(-v^iJm r^(^/^)) = 1 + 0{u-^) 
uniformly for (3 G [0, j3o\. This means that 

4>{v, I3)<iv 

converges uniformly, i.e. given any e > there exists vi such that for any V2> ^i-, 



(A.70) 



4>{v, I3)dv 



< e 



for all P £ [0, /9o]- f{^,P) is differentiable in u, and decreasing, so that 

df 



di 



^ 0. 
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If we let ijj{v,(3) := (j){v' , I3)dv' then integrating by parts gives 



V2 



1^2 



^f{u2,P)e-e 



^2 



dl 



where we have used the mean value theorem for integrals. If additionally, vi > 1 then 
< 1 and we are done. □ 

Corollary A. 4. 

/ / P^(6e"^"''^^^"'"'^d/?de = / / P^(0e"^"''(^^"'"^ded/3 (A.71) 
Jo Jo Jo Jo 

Proof. This follows immediately from lemma lA. 31 see, for example, §11.55.11 of [St]. □ 
Lemma A. 5. The integral 



exp(-v//3r^(0-icT/3)d/3 
is uniformly convergent for G [0, .^o] for all > 0. 

Proof. We, in fact, prove the stronger statement that the integral in question is uniformly 
convergent for ah ^ > 0. Taking M(/3) := e""^*^. 



□ 



|exp(-V/3r^(0-ia/3)| <M(/?) 
and the integral is uniformly convergent by the Weierstrass M-test. 
Lemma A. 6. The iterated integral 



'0 Ji 

converges uniformly for R > 1. 



d/3de 



(A.72) 



Proof. We shall first consider the case where cr < 0. A lengthy calculation gives 



Ji io- id 



-iR^a-RTr,{^) 



2e-3W4(_^)3/2 



exp 



erfc e-'^'/^y/^ 



i7r/4 



erfc Re-'^'/^V^ + 



By Lemma Fa. 11 Trj{S,) ~ as ^ ^ oo, so 

f'OO 

/ P^(Oe-'^'"-^"''(«)de 
Jo 
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is uniformly convergent for i? > 1 by the Wierstrass M-test with M(^) := Ce '^^^^'^ for 
some constant C which does not depend on ^. 

We can write 



exp 



4cT / \ 



= exp {-R^ia - Rt,{0) w (^Re'^/"^ + ^^^^^£^ j 

and since w{z) = 0{z~^) as z ^ oo and |t^(^)| ^ r*, 

w (^Re'-/'^^+^I!^^0^^ = 0(1) 

as ^ — > oo, uniformly for R > 1. Since | exp(— iiT^(,^))| ^ exp(— r^(^)) we see that the 
convergence of 

fpMe.. (^) e* {..e-'"^V^^ ^) 

is uniform for R > 1, by the Wierstrass M-test. 
In the case cr = we have the simpler integral 

exp(-7^r,(e))d/3 (r,(e) (e-^-^) - i?e-^-"«)) 

The integral with respect to ^ then converges uniformly by the Wierstrass M-test, since 

\expi-RT^{0)\ ^ exp(-9le T^(e)) 

and 

|i?exp(-i?r^(0)| ^ ^ exp{-^mt T^iO) 
for R>1. □ 

The following theorem from §11.55.111 of |^ describes criteria which permit the inter- 
change of order of two improper integrals. 

Theorem A. 7. Let f{x,a) be continuous in ai ^ a ^ 02 and c ^ x ^ d, where both 02 
and d may be arbitrarily large, and; 

f{x,a)dx be uniformly convergent for a G [01,02], 

l-OO 

ii) I /(x,a)da be uniformly convergent for x G [c, d], 

poo /■R 

/ / f{x,a)(ladx be uniformly convergent for R [ai, 00]. 

J c J ai 

then 

POO POO POO POO 

/ / f{x,a)dxda= / / f{x,a)dadx. 

J a\ J c J c J ai 
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Applying theorem lA.TI to the integral in (|4.45|) allows us to conclude the following. 
Proposition A. 8. 

/ / P,(Oexp(-v^r,(e)-ia/3)dCd/3= / / P^(0 exp(- v^r,(0 - i(T/3)d/3de. 
JO JO Jo JO 

Proof. This follows from theorem I A. 71 with lemmas IA.31 IA.5I and IA.61 together with 
corollary IA.4I □ 



B Appendix: Simplification of ( 15.551 ) 

We here consider some technical points that arise in section |S1 
Lemma B.l. The integral 



is uniformly convergent for a £ [— i?^,0] for any R> 0. 
Proof. Expanding the w function, using lemma 



2(-a)3/2 \^ j id \T.r^{i? 

as ^ — > OO where the implied constant is independent of cj G [—R, 0]. Since PriiO^i = o 
the leading order term in the expansion of (|B.73|) has zero real part, and the integral of 
the remainder converges since t^(C)~^ ~ (a^)~^ as ^ ^ oo. □ 

Proposition B.2. We have 



Proof. We make the substitution 2p = (— o") to give 



2vu{t)dt 

75, 

where t £ C follows the contour 7^ connecting to — . Since w is an analytic 

function, we can write 



2(-(7)3/2 I 2^^ / \ 2R 
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where W is the antiderivative of w satisfying 

dW 



dz 



w{z) and W{0) = 0. 



By making the substitution ^ i-^ R^, we see that 



33i7r/4 



2R 



2R 



di (B.74) 



where we have, additionally, split the range of integration into two regimes and the first 
integral made the substitution v = ^ . For the first integral in HB.74I) we consider 



R 



exp 



ivr _ lR^\ / e3W%(iV^\ dz^ 
4 ^ ) \ 2i? ) ^ 



which comes from the first term of Prf{R\/v). The second term of Pr^ can be handled in 
the same way. Differentiating, 



d_ 
d^ 



W 



^3i7r/4 



2R 



^3in/4: /p3i7r/4 



4V^ 



W 



2R 



r;(i?^). (B.75) 



dr„ 



Since is bounded for ^ G M, we deduce from ()B.75|) that there exists a constant K 
independent of R such that 



d 

d^ 



W 



' e^'^/\{R^) ^ 



2R 



K 



(B.76) 



Let 



which satisfies 



:= — vvrerfc 
dV' R /ivr iR^v 



We can then use integration by parts, 



ivr \R^u 



R 

2^"nT 4 



w 



' e^'^/\{R^) \ di^ 



2R 



f^l\{R^/u)\ 



2R 







w 



' e^^/\{R^) \ 



(B.77) 

dv 
(B.78) 



as R ^ CO, since 



33i7r/4 



2R 



23 



and 




and the fact that the final integral in ()B.78|) converges uniformly by (|B.76|) . 

For the second integral in ()B.74|) we apply Taylor's theorem and lemma lA. II to get 

This gives 

as -R — > oo. The integral which remains is of a form for which the asymptotic series may be 
derived by the method of repeated integration-by-parts ^BlHaj to see that this contribution 
also vanishes in the limit R ^ oo. □ 
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